Cohesive energies of Fe-based glass-forming alloys 
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We calculate the cohesive energies of Fe-based glass-forming alloys in the B-Fe-Y-Zr quaternary 
system. Our ab-initio calculations fully relax atomic positions and lattice parameters yielding en- 
thalpies of mixing at T=0K. We examine both the known equilibrium and metastable phases as 
well as a selection of plausible structures drawn from related alloy systems. This method, generally 
reproduces experimentally determined phase diagrams while providing additional information about 
energetics of metastable and unstable structures. In particular we can identify crystalline structures 
whose formation competes with the metallic glass. In some cases we identify previously unknown 
structures or observe possible errors in the experimental phase diagrams. 

PACS numbers: 61.50.Lt,61.43.Dq, 71.20.Be, 81.30.Bx 



I. INTRODUCTION 

Calculation of alloy phase diagrams from first- 
principles is necessary to achieve the goal of materials 
by design [|. Bulk metallic glass- forming alloys, which 
often contain three or more chemical elements provide a 
useful test case. Recently discovered many-component 
alloys 0, 0, 3 solidify in amorphous structures at rela- 
tively low cooling rates. These materials display intrigu- 
ing and potentially useful mechanical properties includ- 
ing nearly perfect elasticity @- Amorphous Fe-based al- 
loys are interesting for both their structural and their 
magnetic properties. Achieving bulk glass formation 
could extend the range of potential applications of these 
materials. 

To understand factors limiting bulk glass formation, 
we perform ab-initio total energy calculations on the 
quaternary compound B-Fe-Y-Zr as well as its binary 
and ternary subsystems such as B-Fe and B-Fe-Zr. Our 
calculations use the plane-wave electronic density func- 
tional theory program VASP 0, Q- We identify the 
crystalline structures whose formation competes with 
the amorphous structure of the supercooled liquid. To 
this end, we calculate the cohesive energies of stable, 
metastable and hypothetical crystal structures through- 
out the alloy composition space. Standard metallurgi- 
cal databases [H, ||i [T(|list known stable and metastable 
structures. Chemically similar alloy systems provide hy- 
pothetical structures to test. 

Cohesive energies of stable and metastable phases yield 
thermodynamic driving forces for crystallization. The 
structural complexity of these phases gives some insight 
into the possibility of their nucleation and growth during 
a rapid quench. Our main results are the identification 
of the structure types CFe3 (Pearson symbol oP16) and 
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CeCr23 (Pearson symbol cF116) as the two main com- 
petitors to the B-Fe glass. Alloying with Zr does sur- 
prisingly little to destabilize the CgCr23 structure, while 
alloying instead with Y does reduce stability of this struc- 
ture. On the other hand, alloying with Y stabilizes cer- 
tain other ternary structures. On this basis we deduce 
advantageous composition ranges. 

Energetically favorable structural motifs identified 
within these phases can be compared with structural 
models of the metallic glass. We identify a class of Boron 
coordination polyhedra related to the trigonal prism [Tl| 
with some distortions. These polyhedra may be arranged 
in many ways, some leading to simple crystal structures 
with very low energy but many more whose energies and 
local structure closely resemble the glass. We call the en- 
tire class of these structures "amorphous approximants" , 
by analogy with the concept of approximants to qua- 
sicrystal structures [T3 | . 

The following two sections of this paper (sections [D] 
and present our calculational methods and the re- 
sulting cohesive energy data. Although we study here Fe- 
based glass-forming alloy systems, our basic method can 
be applied to any alloy system. Besides checking known 
experimental phase diagrams, our calculations provide 
energetic information which is often not known experi- 
mentally, especially in the case of metastable and amor- 
phous structures. In addition, we propose likely struc- 
tures for compounds whose existence was known but 
whose structures were unknown, for example, B4Fe4Y, 
B4FeY and BgFe2Y5. Conversely, in some cases our re- 
sults call into question details of the established phase 
diagrams. For example: the claimed high temperature 
stability of BZr is most likely only metastability in real- 
ity; the C07Y2 structure may be stable in Fe-Y though 
it has not been reported; the phase FeiyYo, related to 
important permanent magnet materials [l3| . is possibly 
only a high temperature phase. Further, we can predict 
phase diagrams of alloy systems such as B-Y-Zr, Fe-Y-Zr 
and B-Fe-Y-Zr that have not been established experimen- 
tally. 
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Finally, section llVl analyzes the crystal structures and 
correlates their atomic arrangements with their cohesive 
energies. Amorphous approximants are presented in sub- 
section |IVD] 



II. METHODS 

Our interest in the binary B-Fe and ternary B-Fe-Zr 
and B-Fe-Y compounds led to the study of all elemen- 
tal, binary, ternary and quaternary combinations of the 
elements B-Fe-Y- Zr. We selected structures for study 
that are known as stable or metastable structures in the 
phase diagrams of these alloy systems or chemically sim- 
ilar alloy systems. For example we consider known C-Fe 
structures (e.g. CFe3.oP16) with B replacing C, etc. Our 
notation for structure type is to first give the prototype 
(some familiar isostructural compound, e.g. CFes) fol- 
lowed by the Pearson symbol (indicating the point sym- 
metry, centering information, and number of atomic sites 
per unit cell, e.g. 0PI6 for orthorhombic Primitive 16- 
atom cell). 

Our sources for established phase diagrams and struc- 
tures include standard references 

B la E3, 

individual 

publications and private communication. Information 
from these sources has been compiled into a database 
containing over 1000 structures that we search to match 
criteria such as chemical elements, stoichiometry, and 
atomic size ratios. Some additional structures examined 
are liquid and amorphous structures, obtained from ab- 
initio molecular dynamics simulation. 

Our ab-initio calculations use the program VASP (ver- 
sion 4.5.5) together with the projector-augmented wave 
method, an all-electron generalization of the pseudopo- 
tential approach 0, 0|. We employ the Perdew- 
Wang generalized gradient approximation 0] (GGA) 
exchange-correlation functional with the Vosko- Wilkes 
Nussair s P m interpolation. These choices give ex- 
cellent results for bulk elemental Fe [15J . GGA is needed 
instead of LDA is necessary to properly reproduce mag- 
netization and lattice constants [18| . Our magnetic cal- 
culations are spin-polarized (i.e. collinear magnetization) 
and are employed for any structure containing 50% Fe or 
higher. 

VASP solves for the self-consistent electronic structure 
in reciprocal space, using a plane wave basis. It requires 
that we choose the reciprocal space grid appropriately 
and demonstrate convergence in the number of k-points 
used and the plane wave energy cutoff. We construct k- 
point grids whose spacing is nearly isotropic in reciprocal 
space. Mostly we use Monkhorst-Pack grids, although for 
hexagonal structures we use Gamma-centered grids. Our 
k-point density is sufficient that all structural energies are 
converged to a precision of 10 meV/atom or better. All 
energies for structures that lie on or near the convex hull 
are converged to a precision of 1 meV/atom or better. 

In general our relaxations allow variations of cell vol- 
ume and shape, as well as atomic displacements, consis- 



tent with the symmetry of the starting structure. Re- 
laxations run until an accuracy of lmeV/atom or better 
is reached. During relaxation we use Methfessel-Paxton 
Fermi-surface smearing with width 0.2 eV (the VASP 
default choice [l^. When smearing is employed we re- 
port the energy (extrapolated to zero smearing), not the 
fictitious free energy. For many structures on or near 
the convex hull, we recalculated the energy of our best- 
relaxed structure using the tetrahedron method without 
smearing. This test confirms we reached our 1 meV/atom 
precision goal. 

Table H] demonstrates convergence with respect to k- 
point grid and plane wave cutoff energy, reporting the 
energies of two metastable variants of Fe 3 B, one with 
structure type Fe3C.oP16 and one with structure type 
Ni3P.tI32. We see how energies E and energy differences 
AE converge as the k-point mesh grows at fixed energy 
cutoff. The energy difference AE at low precision has 
even the incorrect sign, while medium and high precision 
agree to within the desired lmeV/atom accuracy. For 
the results presented below we employ a constant energy 
cutoff of 319 eV, consistent with medium precision. 

Certain calculated structural quantities can be com- 
pared directly with experiment. For example, for 
Fe3B.oP16 we obtain converged volume 9.78 A 3 /atom, 
b/a = 1.23 and c/a = 1.52 compared with the experimen- 
tal values, respectively, of 10.09 A 3 /atom, 1.22 and 1.51. 
Likewise, for Fe3B.tI32 we find volume 9.61 A 3 /atom 
and c/a = 2.03 compared with experimental values 
10.06 A 3 /atom and 2.01. Our underestimation of the vol- 
ume reflects both thermal expansion (experimental vol- 
umes are at room temperature, while our calculations are 
for T=0K ground states) and known systematic errors 
associated with density functional theory. 

The composition space of an V-component alloy is a 
set of N composition variables {xi : i = 1,2, N} obey- 
ing y^^L 1 Xj = 1. The set forms an N — 1 dimensional 
simplex (respectively, a point, line segment, triangle and 
tetrahedron for N — 1, 2, 3, 4). Structural energies form 
a scatter-plot over this simplex. Stable low temperature 
phases lie on vertices of the convex hull of the energy 
versus composition scatter-plot. Edges and facets of the 
convex hull represent coexistence regions of the phases at 
adjoining vertices. Lines and triangles joining low tem- 
perature phases phases will be referred to as tie-lines and 
tie-planes. A tie-surface in general refers to the hyper- 
plane joining N or fewer points in the V-component en- 
ergy scatter-plot. 

The tie-surface connecting all pure elements in their 
lowest energy structures forms a useful reference for al- 
loy energies. The distance AHf or of an alloy energy from 
the tie-surface joining pure elements is known as its en- 
thalpy of formation (enthalpy because volume relaxation 
means we work at fixed pressure, P = 0). Strong com- 
pound formation is reflected in large negative enthalpy 
of formation. 

High temperature phases should lie above the convex 
hull, but be sufficiently close that entropic effects (e.g. 
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TABLE I: Convergence of magnetization and energy with re- 
spect to k-point density and plane wave cutoff. Numbers in 
left-hand column are the product of lattice parameter times 
the number of k-point subdivisions. Values of E c are plane 
wave cutoff energies in eV 
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phonons or chemical substitution) can stabilize them. 
Metastable phases also should lie close to the convex hull, 
so that their free energy is less than the liquid free en- 
ergy at temperatures below freezing. Although AHf or 
is usually negative for high temperature and metastable 
phases, their energy difference AE from the convex hull 
is small and positive. The value of AE is a measure of 
the thermodynamic driving force for decomposition into 
the appropriate combination of stable phases. 

Using these methods, we built a database of structural 
energies. For a given A^-component alloy system of inter- 
est we extract from our database energies of structures 
containing all, of just some, of the chosen elements. We 
use a standard convex hull program (qhull [IS]) to iden- 
tify stable structures and the coexistence regions that 
connect them. Based on the output of this program, we 
calculate values of AHf or and AE for every structure. 

Our methods introduce systematic errors associated 
with the PAW implementation and even with the un- 
derlying density functional theory. Provided these errors 
vary smoothly with composition, the identity of convex 
hull vertices will not be affected in most cases. However, 
the tie-lines, tie-planes, etc. grow progressively more sen- 
sitive to error. It is probable that even when we correctly 
identify the stable phases, we may misidentify their co- 
existence regions. 

Previous workers carried out analogous studies, though 
mainly on binary alloys. Hafner |2l| wrote a gen- 
eral introduction to the subject of ab-initio alloy phase 
diagram prediction. Miedema, de Boer and cowork- 



ers [22j performed extensive semi-empirical analysis of 
binary alloy systems, including almost all binary alloys 
of Fe ■ Others 0, |2^, |2|| take more rigorous ap- 
proaches, sometimes including finite temperature effects 
of vibrational and configurational entropy. Many co- 
hesive energy calculations have been collected in online 
databases [23, |2il • All the data we present here, and 
a great deal more, can be found on the WWW |30j . 

III. RESULTS 

The N = 4-component B-Fe-Y-Zr quaternary alloy 
system contains many subsystems: four pure elements, 
six binary alloy systems and four ternaries. This sec- 
tion presents our results in order of increasing number 
of components. We adopt alphabetical order in naming 
alloy systems, because this brings some order to the pro- 
liferation of chemical combinations in multicomponent 
systems. 

In most cases we reproduce the known equilibrium and 
metastable phase diagrams with surprising accuracy. In 
a few cases lingering discrepancies may reveal limitations 
of our method or cast doubt on the accepted phase di- 
agrams. In some cases where the existence of a phase 
was known experimentally but not its structure we sug- 
gest probable structures. Our quaternary and two of our 
ternary phase diagrams have never been determined ex- 
perimentally. 

A. Pure elements 

Each of the four elements under study exhibits solid- 
solid phase transitions in addition to its melting transi- 
tion, and verifying the relative energies of the different 
structures is a nontrivial test of our calculational method. 

The precise low temperature structure of Boron is un- 
known. The presumed equilibrium phase, designated f3 
is rhombohedral with approximately 108 atoms per unit 
cell. Atomic positions are known, but there is partial 
occupancy and probably strong correlations among the 
partially occupied sites, that have not been adequately 
resolved. 

The basic structure of /3-Boron (denoted B.hR105) 
consists of overlapping B156 clusters with icosahedral 
symmetry, located at the vertices of the primitive rhom- 
bohedral cell [33. The 12 * (12 + 1) = 156-atom clus- 
ter assembles twelve 12-atom icosahedra surrounding one 
central icosahedron. All icosahedra all empty. A single 
extra atom at the body center of the rhombohedral cell 
is not part of any B156 cluster. An alternate description 
concentrates on non-overlapping Bs4 clusters, obtained 
by removing the outer halves of the 12 outer icosahedra 
of B156. In addition to the B§4 cluster the basic structure 
contains two Bio clusters connected by the extra B atom 
in the center of the rhombohedral primitive cell, yielding 
105 atoms per unit cell. 
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Wyckoff positions of B.hR105 are labeled B1-B15, 
with the B15 site the cell center. Structural refine- 
ments ji^l 0s2, 03 0] fi 11 ^ t ne B13 sites surrounding the 
B15 atom at the center of the rhombohedral cell only 
fractionally occupied, while an extra B16 Boron atom 
approximately compensates the missing electron density. 
The refinement by Slack |34j additionally reports five 
other Boron sites with small occupancy factors, and pro- 
poses tentative a model for occupancy correlations. The 
Pearson symbol for the Slack model is hR141. 

Our calculations show the basic B.hR105 structure 
higher in energy than the presumed metastable a-B 
structure (B.hR12) by 25 meV/atom. We find, in agree- 
ment with Slack [34( but contrary to naive assumption, 
that the choice of which B16 atom to insert is not corre- 
lated to the B13 site occupancy. The B16 atom lowers the 
total energy by as much as 8 eV per primitive rhombohe- 
dral cell, narrowing the energy difference between a and 
(5 boron to ~11 meV/atom. Insertion of a second B16 
atom, or other more complex modifications, are second- 
order corrections to the total energy. Our best model of 
the fractional occupancies additionally replaces one B13 
atom by a new B19 atom (106 atoms per primitive cell), 
provided the B19 atom is not a nearest neighbor of the 
B16 atom. This model remains unstable relative to a-B 
by ~ 5 meV/atom. 

Pure elemental Iron passes through four solid phases, 
a-8, as temperature rises from T=0K to melting, a is 
BCC (Pearson cI2) and ferromagnetic. At its Curie tem- 
perature it transforms to (3 which is also BCC. At higher 
temperatures it transforms to 7, which is FCC (Pear- 
son cF4) and finally to 8, again BCC, before melting. 
First-principles calculations 35] show that the magnetic 
ground state of 7-Fe is a noncollinear antiferromagnet, 
while at the experimental atomic volume it is a collinear 
antiferromagnet |36j| . According to our calculations, the 
relaxed energy of this collinear antiferromagnetic struc- 
ture is 81 meV/atom above the energy of a-Fe. 

The low temperature a phases of Yttrium and Zirco- 
nium are both hexagonal (hP2) and their high tempera- 
ture (3 phases are both BCC (cI2). Our calculations agree 
with these facts. 



B. Binaries 

1. Fe-B 

The established Fe-B phase diagram contains just two 
compounds, Fe2B.tI12 and BFe.oP8, each with sim- 
ple structures. A number of metastable phases exist 
in the Fe-rich end, namely CFe3.oP16, Ni3P.tI32 and 
CeCr23.cF116. These occur close to the deep eutectic 
at 17 % B, and hence are important competitors to glass 
formation. 

Our calculations shown in Fig. ^ reproduce the estab- 
lished phase diagram perfectly, with the known stable 
phases lying on the convex hull and the known metastable 



phases lying within 20 meV/atom above (at the eutec- 
tic melting temperature of T=1447 K, a characteristic 
thermal energy is ksT — 125 meV/atom). Structures 
occurring as stable or metastable phases in other alloy 
systems (but not in Fe-B) lie further above the convex 
hull, except for the BRe3.oC16 and C0SC3.0P32 struc- 
tures. On the basis of our calculation, BRe3 and C0SC3 
should occur as metastable phases in Fe-B, although they 
have not been reported. 

A prior ab-initio study of the Fe-B structures investi- 
gating magnetism and bonding [37j reported a covalent 
character to Fe-B and B-B bonds while the Fe-Fe bonds 
are metallic in nature, together with a significant charge 
transfer from Fe to B. Spin polarization of the B atoms 
is weak and opposite to the Fe atoms. Our calculations 
support those conclusions. 

Given our perfect agreement with the established 
phase diagram, it is surprising that our calculated en- 
thalpies (respectively -368 and -308 meV/atom for FeB 
and Fe2B) differ greatly from measured values (respec- 
tively -676 and -707 meV/atom for FeB and Fe 2 B |3S^. 
This may be due in part to the fact that our calcula- 
tion are performed at T=0K while the measurements 
were done at T=1385K. Pure elemental Fe undergoes two 
phase transitions (one structural, one magnetic) as tem- 
perature drops, which could be partly responsible for this 
discrepancy. 

One hypothesis on the glass-forming ability of Fe- 
B is that very simple, easy to nucleate and grow 
crystal structures, are destabilized by the size con- 
trast |13 of Fe (nominal diameter=2.48A) and B (nom- 
inal diameter=1.80A). We observe this principle in ac- 
tion in the Fe-rich end of this phase diagram. Consider 
a substitutional solid solution of Fe and B, around com- 
position BFe3, based on the BCC structure of Fe. The 
Fe3Si.cF16 structure is a particular realization of such 
a structure, in which the Fe and B atoms arrange at 
regular positions. However, the energy of cF16 is much 
higher than the metastable 0PI6 structure, which can be 
reached through distortion of the cF16 lattice. It seems 
that the lattice strain caused by size mismatch destabi- 
lizes the BCC solid solution, converting cF16 into 0PI6. 

We note in addition that B will not stabilize the FCC 
structure of Fe by substitution to form the AuCu3.cP4 
structure, though such stabilization does occur with 
larger atoms, for example FesGe. 

Alternatively, B might enter as an interstitial, as 
indeed C enters into FCC lattices of Fe in octahe- 
dral or tetrahedral sites (respectively, in the metastable 
Fe3C.hP8 and Fe4C.cP5 structures). However, owing to 
the slightly larger size of B compared with C (nominal 
diameter=1.43A), these structures are far above the con- 
vex hull in the B-Fe energy scatter plot. 

It seems that only fairly complicated crystal structures 
exist near compositions of about 25% B. The difficulty of 
nucleating and growing these structures may aid in glass 
formability. An estimate of the thermodynamic driving 
force for nucleation can be obtained by comparison of 
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FIG. 1: Enthalpies of formation and their convex hulls for the B-Fe binary alloys. Notation: heavy circles denote known low 
temperature phases, light circles denote known high temperature phases, diamonds denote known metastable phases, squares 
denote unreported structures. 



three energy scales. The metastable structures are about 
10-20 meV/atom above the tie-line joining Fe2B to a- 
Fe. Further details of Fe-B structures, especially focusing 
on B-atom environments and the occurrence of trigonal 
prism structures is given in the discussion section llVl 

A quenched amorphous structure at composition 
FesoB 2 o is about 170 meV/atom above the tie-line, and 
the liquid is 350 meV/atom above. Candidate liquid 
and amorphous structures were produced by liquid state 
molecular dynamics on 100 atom samples using Nose dy- 
namics at temperature T=1500K, then quenched by con- 
jugate gradient relaxation. Atomic displacements under 
relaxation are large, averaging about 0.6 A. The liquid 
runs were at fixed volume, 7% greater than the volume of 
the quenched amorphous samples. All molecular dynam- 
ics and quenching runs were done with spin polarization, 
using the T k-point only. 



2. Fe-Zr and Fe- Y binaries 

A promising glass-forming strategy is to start with a 
good glass-forming binary such as B-Fe, then add one 
or more additional elements to further destabilize any 
crystalline structures. Large atoms such as Zr (nominal 
diameter=3.18A) and Y (nominal diameter=3.55A) are 
promising because they associate well with Fe and B but 
they differ strongly in size from either Fe or B. Before 
turning to these ternary and quaternary compounds, we 



briefly examine the Fe-Zr and Fe-Y binaries. 

The established Fe-Zr phase diagram |4(j con- 
tains three low temperature compounds, Fe2Zr.cF24, 
FcZr2.tI12 and FeZr3.oC16, one high temperature com- 
pound Fe23Zr6.cF116 (this cF116 structure is quite dis- 
tinct from the C6Cr23.cF116 structure which confusingly 
shares the same stoichiometry and Pearson symbol) and 
one metastable compound Fc 2 Zr.hP24. Our calculation 
(see Fig. |5J) is in excellent agreement with experiment. 
Every known low temperature phase lies on the con- 
vex hull, and the high temperature and metastable com- 
pounds lie close above it. 

One unknown structure, Fe3Y.hR12, appears on the 
convex hull where no stable compound is known exper- 
imentally. Most likely this reflects an inaccuracy of our 
methods. Because its stability relative to the tie-line join- 
ing Fe2Zr to pure Fe is about 1 meV/atom, small errors 
(either calculational or arising from approximations of 
density functional theory) could account for this differ- 
ence. Alternatively, the phase could truly be stable, but 
hard to observe experimentally because the driving force 
for its formation is weak. This matter requires further 
theoretical and experimental analysis, but for the study 
of glass formation all we care about is that its energy lies 
close to the tie-line from Fe2Zr to a-Fe. 

Now consider Fe-Y. The established phase diagram 
contains all the same phases as Fe-Zr and two addi- 
tional phases, Fe3Y.hR12 and FeiyY2. Certain features 
of the phase diagram are thermodynamically improb- 
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able |4ll l42j . The close proximity of line compound 
Fe23Yg.cF116 (at low temperature) to the line com- 
pound Fe 3 Y.hR12 is highly unlikely. We presume that 
Fe-23Y 6 .cF116 is stable at high temperatures only, where 
the phase diagram shows a broad composition range. The 



strong asymmetry of the liquidus of Fe3Y.hR12 also is im- 
probable but we have no proposed alternative at present. 

The structure of Fei7Y 2 has not been properly iden- 
tified, and is believed to occur in at least two vari- 
ants. We follow the lead of Massalski and Okhamoto 
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and identify these variants as Th2Zni7.hR19 (low tem- 
perature) and Nii7Th2.hP38 (high temperature). Other 
reported variants [8j of this phase are Fei7Ho2-hP44 and 
Ni 19 Th 2 .hP80. The occurrence of several structural vari- 
ants, most with partial occupancy, suggest a possible 



entropic stabilization mechanism by structural disorder. 
See the synopsis of the basic structure and its degrees of 
freedom in the discussion section llVl 

Our calculations for Fe-Y (Fig. |3J) present certain dis- 
agreements in comparison with experiment. Notably, all 
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variants of Fei7Y2 lie above the convex hull and thus are 
predicted as high temperature or metastable. Meanwhile, 
Fe5Y.hP6, believed metastable, touches the convex hull 
as does the unreported structure C07Y2.hR.l8. 

For the failure of Fei7Y2 to meet the convex hull, 
three possible explanations are: (1) The structure of 
Fei7Y2 has not yet been correctly determined (we men- 
tioned this above and discuss it further in section TlV Bj l: 
(2) Our calculations are seriously flawed and unable to 
properly compare the energies of Fe-Y compounds (we 
checked that changes in cutoff energy, pseudopotcntials 
and exchange-correlation functional have no significant 
impact); (3) Fei7Y2 is only metastable or high tempera- 
ture and the true low temperature state is a coexistence 
of FesY and pure Fe. 

A simple mechanism to explain high temperature sta- 
bility is note that the energy of Fei7Y 2 lies below the 
tie-line from FesY to 7-Fe. Because 7-Fe is the phase 
with which Fei7Y2 coexists from melting down to about 
T=900C, it may be difficult to observe decomposition of 
Fei7Y2 at low temperatures where it competes instead 
with a-Fe. This scenario suggests the possibility that 
other Fe-based alloy phase diagrams could be incorrect 
at low temperatures, which could have significance for 
the engineering of magnetic materials [l3|. 

The ab-initio calculations of total energies for the Fe- 
Zr and Fe-Y alloys were difficult because of the compli- 
cated magnetic properties of the Fe-rich structures. We 
mention two important observations. (1) The Zr and Y 
atoms have magnetic moments pointing opposite to the 
Fe atoms. Moments are typically in the range of +1.8 to 



+2.4 hb for Fe and in the range of -0.2 to -0.4 hb for 
Y or Zr in spin polarized calculations. Similar values for 
Fe-Y alloys were reported in prior calculations j43| and 
experiments 0|. (2) Magnetism couples strongly with 
atomic volume leading to multiple self-consistent solu- 
tions of the DFT. Generally one finds: a nonmagnetic, 
low volume, high energy solution; a strongly magnetic, 
high volume, low energy solution; occasional additional 
solutions of intermediate magnetism, volume and energy. 
Presumably this is related to the strong magnetovolume 
effects that actually occur in Fe-rich compounds 0, \^ . 

Comparing our calculated enthalpies of formation with 
published experimental data |3Sj , we again find that our 
calculated values lie well below the published data, in 
the case of Fe-rich alloys, probably as a result of the high 
temperatures at which the experiments were carried out. 
At lower Fe-content our data is fairly consistent with the 
experimental data. 



3. Other binaries: B-Zr, B- Y and Zr- Y 

Next we turn to the B-Zr and B-Y binaries. The es- 
tablished B-Zr phase diagram exhibits three compounds. 
The well known stable phase B2Zr.hP3 is very strongly 
bound. The other two phases exist only at high tem- 
peratures: Bi2Zr.cF52 melts congruently, while BZr.cF8 
exists only over an intermediate temperature range be- 
low all liquidus temperatures 0, |4g • Our calculation 
(Fig.^J supports stability of B2Zr and high temperature 
stability of B^Zr, but strongly contradicts the existence 
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of BZr.cF8. Indeed, at this composition the best struc- 
ture we find is BCr.oC8, but that too is highly unsta- 
ble. We explored nearby compositions, and even added 
traces of C (CZr.cF8 is a stable compound in the C- 
Zr binary system), but we cannot find any structure 
within a reasonable distance of the convex hull. Prob- 
ably BZr.cF8 is a metastable structure formed during 
rapid quench |4*5Ll46j. 

Comparing our calculated enthalpies of formation with 
published data [3^| we find excellent agreement. For 
B12Y we calculate -219 meV/atom compared with the 
experimental value -213 meV/atom. The experimental 
data is taken at a fairly low T=298K. For B2Zr we calcu- 
late Aif/ or =-999 meV/atom compared with the experi- 
mental value -1074 meV/atom. 

The enthalpy of B2Z1- reflects strong covalent B-B 
bonding. Densities of states among transition-metal 
diboridesji]} exhibit a strong pseudogap associated with 
Boron p-states. As one moves across the transition-metal 
series, the Fermi level falls in the pseudogap for group IVa 
elements (Ti/Zr,Hf), leading to strong peaks in cohesive 
energy. This is also the likely cause of the wide compo- 
sition range of B 2 Nb - since Nb lies just to the right of 
Zr and the Fermi energy of i? 2 Nb lies just to the right of 
the pseudogap, Nb vacancies can move the Fermi energy 
closer to the gap, resulting in a low vacancy formation 
energy. 

In contrast to B-Zr, B-Y contains several B-rich 
phases, including B2Y.I1P3 and Bi2Y.cF52 as in B-Zr, 
and additionally B4Y.tP20 and BgY.cP7. The last one is 
of uncertain composition 9], with the experimental dia- 
gram showing a composition range at low temperatures, 
contradicting the notion that alloys should reach defi- 
nite compositions as T — > OK. We investigated Yttrium 
vacancies within a 2 x 2 x 2 supercell of BgY.cP7 and 
found that removal of a single Y out of 8 was favorable 
and lowered the value of AE to +15 meV/atom. Finally, 
there is a phase Bg6Y.cF1880 of whose gigantic unit cell 
size prevents us from calculating cohesive energy. Our 
calculated convex hull (Fig. |5J) agrees perfectly with the 
experimental data except in the case of B@Y which we 
find is unstable at low temperature. Owing to the lack 
of a definite low temperature stoichiometry, B$Y most 
likely is unstable at low temperatures [41], |43 ■ No exper- 
imental data is available for enthalpies of formation. 

Finally, consider the Y-Zr binary (Fig.[f)J). The estab- 
lished phase diagram contains no compounds. We inves- 
tigated several possibilities appropriate for their atomic 
size ratio and found none stable. 



C. Ternaries 

1, B-Fe-Zr and B-Fe-Y 

The ternary alloy system B-Fe-Zr exhibits no known 
ternary compounds. Our calculations (Fig. 01 generally 
support this, but we do identify one stable compound, 



with structure, B4CrY.oP24. We convert the binary 
structure C6Cr 2 3.cF116 to a ternary by substituting the 
large Zr atoms on sites that have the largest Voronoi vol- 
ume. These turn out to sites of Wyckoff type 8c, resulting 
in the intrinsic ternary structure type BgCo2iZr2.cF116. 
Interest in this structure is motivated by its proximity 
to the glass-forming composition, and the metastability 
of this structure is discussed later in section IIVI The 
quenched structure is, as before, a 100 atom model metal- 
lic glass reached by molecular dynamics and quenching. 

The reason that B-Fe-Zr exhibits so few (i.e. just one) 
stable ternary compounds is that the enthalpy of for- 
mation of B 2 Zr is very large, as discussed above in sec- 
tion Covalent bonding of B2Zr is so strong that 
even Fe-rich systems find it advantageous to phase sep- 
arate into a mixture of B2Zr plus Fe alloyed with which 
ever of Zr or B remains in excess. 

In contrast, the B-Fe-Y ternary alloy system (Fig. [SJ 
exhibits many stable ternary compounds, because the 
bonding of B2Y is less strong than B2Zr. Those com- 
pounds with known structures are the Fe-rich compound 
BFei 4 Y 2 .tP68, and in the Fe-poor region B 2 Fe 2 Y.tI10, 
B 6 Fe3Y4.hR13, ByFeY 3 .oC44 and a metastable struc- 
ture Bi4Fe62Y3.cI158. All previously known stable 
structures touch the convex hull. Additionally there 
are stable compounds of unknown or partially known 
structure at B3FeY2, B4Fe4Y and B4FeY. Our calcu- 
lations reveal these compounds to take the structure 
types, respectively, B 6 Fe 2 Nd5.hR13, B 4 Co4Nd.tP18 and 
B 4 CrY.oP24. 

Surprisingly, we find a previously unknown compound 
of structure type BCeCo4.hP12 on the convex hull in 
the Fe-rich region. Nearby we find structure type 
B2Nd3Nii3.hP18 just slightly higher in energy. It would 
be of interest to explore these compositions experimen- 
tally in more detail. 

At ^80% B content, we find three stable crystals in 
the B-Fe-Y ternary: B 4 CrY.oP24, B 6 ReY 2 .oP36 and 
B 7 ReY 3 .oC44. One of these, B 4 CrY.oP24, is also sta- 
ble in the B-Fe-Zr ternary. All three of these structure 
types can be considered as approximants to decagonal 
quasicrystals. However, we have not identified systematic 
extensions toward truly quasiperiodic structures, and we 
are not prepared to predict the occurrence of decagonal 
quasicrystals in these compounds. At present no B-based 
quasicrystals are known. Further discussion can be found 
in Ref. |4S| . 

One structure reported in the B-Fe-Y system, 
Bi 4 Fe 62 Y 3 .cI158, has a very high energy (AE = 330 
meV/atom) and large initial forces (as high as 1.7 eV/A). 
Even after large atomic displacements during relaxation 
the energy remains very high. We believe the experimen- 
tally reported structure is incorrect. 

We are impressed by the faithfulness with which our 
calculations reproduce systematic differences in the phase 
diagrams of Y- and Zr-containing alloys. Despite their 
adjacency in the periodic table, and the consequent sim- 
ilarities in atomic size, electronegativity and preferred 
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FIG. 7: Convex hull and metastable phases of the B-Fe-Zr ternary system. For binary structure types see Figs. 111401 Plotting 
symbols as in Fig. Energy units are meV/atom. 



structure types, those details on which the accepted 
phase diagrams do differ are almost always correctly re- 
produced. 



8. B-Y-Zr and Fe-Y-Zr 

The B-Y-Zr ternary diagram has not been experimen- 
tally determined. We have explored it using the methods 
described above. The only stable ternary compounds we 
find (Fig. I1U[) are extensions of certain binaries into the 
ternary. Notably B2(Zr,Y).hP3 exhibits complete misci- 
bility of Zr and Y in this pseudobinary structure. Addi- 
tionally Bi2Y.cF52 extends part way into the ternary. 

The Fe-Y-Zr ternary diagram has not been experimen- 
tally determined either. Our calculations (Fig. |3J) sug- 
gest that Fe 2 (Y,Zr).cF24 and Fe 3 (Y,Zr).hR12 both ex- 
tend across the full ternary diagram, but no other bina- 
ries appear to extend far into the ternary. 



D. B-Fe-Y-Zr Quaternary 

The B-Fe-Y-Zr quaternary has not been experimen- 
tally determined. No quaternary structures are reported 
in the standard references. We have calculated enthalpies 
of formation for 15 different compounds (11 structure 
types, some with alternate chemical occupancies) and 
find no stable quandaries. Our lowest energy structures 
are listed in Table [H] The nearest we come to stability is 
for the structure type B4CrY.oP24, for which AHf or is 
around 3-4 meV/atom for all substitutions of Y and Zr. 
Thus it is likely that the entire Y/Zr substitution yields 
equilibrium structures at high temperatures. We find 4 
meV/atom for substitution of Zr for one of the two Y 
in BFei4Y2.tP68, suggesting significant Zr solubility at 
high temperatures. 

Several factors contribute to the lack of stable B-Fe- 
Y-Zr quandaries: incompatibility of Y and Zr atoms (see 
Fig. ISJ destabilizes quaternaries that are rich in Y or 
Zr; the difficulty of accommodating the slightly differ- 
ing atomic sizes into the same crystal lattice site classes 
destabilizes quaternaries in which Y and Zr are minority 
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3 B 6 ReY 2 .oP36 (-604) 
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21 BFe 14 Y 2 .tP68 (-115) 



unstable (AE>0) 



1 B 4 MoTh.oC24 (2) 

5 B 6 Er 2 Ni 3 .oC22 (89) 

6 B 4 Ir 3 Zr.hP16 (266) 
7B 6 Fe 3 Y 4 .hR13 (150) 
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14 quench (186) 
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FIG. 8: Convex hull and metastable phases of the B-Fe-Y ternary system. For binary structure types see Figs. lll5l3l Plotting 
symbols as in Fig. Energy units are meV/atom. 



TABLE II: Quaternary data 



B 


Fc 


Y 


Zr 


AE 


AH f or 


Structure 


Comments 


67 


17 


4 


12 


3.1 


-653 


B 4 CrY.oP24 


3Zr 


67 


17 


8 


8 


4.4 


-636 


B 4 CrY.oP24 


3Zr 


67 


17 


12 


4 


3.8 


-620 


B 4 CrY.oP24 


lZr 


7 


82 


10 


1 


5.2 


-117 


BFei 4 Nd 2 .tP68 


lZr 


7 


82 


6 


6 


21.9 


-120 


BFei 4 Nd 2 .tP68 


4Zr 


64 
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23 
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-623 


B 7 FeY 3 .oC44 


Zr on 4c 
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23 
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28.4 


-620 


B 7 FeY 3 .oC44 


Zr on 8f 


17 


67 


8 


8 


23.2 


-264 


BCeCo 4 .hP12 




21 


72 


3 


3 


45.9 


-210 


C 6 Cr 2 3.cF116 


Y/Zr on 8c 


20 


70 


5 


5 


211.9 


-61 


quench 





species; the very strong bonding of B with Zr destabilizes 
quaternaries that are B-Zr-rich. 



IV. DISCUSSION 

A. Boron atom environments 

1. Trigonal prisms in Fe-rich B-Fe stable and metastable 
systems 

Trigonal prisms ^l| place six large (Fe) atoms at their 
six vertices, their rectangular (nearly square) faces are 
capped by an additional three large atoms, and they are 
centered by a small (B) atom. They are well known struc- 
tural motifs in compounds with significant contrast in 
atomic size, in the large-atom-rich composition range. 
Too large a size contrast is unfavorable, as they do not 
occur in B-Y or B-Zr binaries. A stringent definition of 
the trigonal prism uses the (radical-planes) Voronoi con- 
struction. The Voronoi polyhedron of the central atom 
should have no triangular face, 3 rectangular and 6 pen- 
tagonal faces. This polyhedron is denoted (0,3,6) in the 
(713,714,125...) notation of Watson and Bennett 49] . 

One stable B-Fe compound, BFe.oP8 which we find 
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5Fe 3 Y.hR12(-192) 
12 Fe 2 Zr.cF24 (-290) 



unstable (AE>25) 
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7 Fe 4 YZr.cF24 (12) 
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10Fe 2 Zr.hP24 (14) 
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1 Th 2 Zn 17 .hR19 (43) 

2 Fe 5 Y.hP6 (36) 
3Fe 23 Zr 6 .cF116(68 
6Fe 2 Zr.hP24(171) 

13 FeZr 2 .tI12 (183) 
14FeZr 3 .oC16(99) 
15 FeZr 3 .oC16 (34) 
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FIG. 9: Convex hull and metastable phases of the Fe-Y-Zr ternary system. For binary structure types see Figs. I2I31H1 Plotting 
symbols as in Fig. Energy units are meV/atom. 



marginally more stable than the BCr.oC8 prototype, 
contains trigonal prisms. In both structures (and also 
in closely related CaCu.mP20 and CaCu.oP40 structure 
types) Fe prisms share 2 out of 3 rectangular faces with 
neighboring prisms, while the third rectangular face is 
capped by an Fe atom. The prisms form columns along 
the shortest periodic direction. The structure has re- 
spectable packing fraction (greater than 0.73) when Fe/B 
atoms are replaced by hard spheres with radius ratio 1.55, 
optimizing the packing fraction. 

In the Fe-rich portion of the B-Fe system (see sec- 
tion 1111 B P , we find seven Fe-rich structures that are 
unstable by less than 50 meV/atom relative to the con- 
vex hull. These include the known metastable phases 
BFe 3 .oP16, BFe 3 .tI32 and B 6 Fe 23 .cF116, and other 
structures BRe 3 .oC16, C 2 Fe 5 .mC28, C 3 Fe 7 .oP40 and 
C 3 Fe 7 .hP20. Of these, the BRe 3 .oC16 and CoSc 3 .oP32 
structures are sufficiently low in energy and differ suffi- 
ciently in composition from the nearest stable crystalline 
phases that we expect they could also occur as metastable 
phases. At higher B-content we find B 3 Ni4.oP28 at low 
energy and possibly metastable. 



With the exception of the cF116 structure, all B 
atom environments in the above-mentioned structures 
are proper trigonal prisms. In the oP16 and oC16 struc- 
tures all Fe atoms are structurally similar, each with 3 
B neighbors. The oC16 structure is characterized by a 
unique stacking mode of the trigonal prisms, forming un- 
terminated columns along the shortest-period a-axis and 
sharing triangular faces. Each Fe has 2 B neighbors, 
and each Fe is simultaneously the vertex of one prism 
and a capping atom of another prism, shifted by a/2. 
The B 3 Ni4.oP28 structure combines the same building 
blocks found in BRe 3 .oC16 and BFe.oP8 - columns of 
trigonal prisms stacked into columns either sharing tri- 
angular faces (oC16) or rectangular faces (oP8). 

The stable phase BFe2 , like a number of other B-TM2 
binary systems, crystallizes into the Al2Cu.tI12 structure 
type. Viewed parallel to the shortest (c) axis, the struc- 
ture is built by two flat layers of Fe atoms each forming a 
nearly regular square-triangle tiling pattern |4S[ . Boron 
atoms occupy interstitial octahedral sites in the network. 
This topology is not optimal for packing atoms with very 
different sizes, so the atomic sizes must not play an im- 
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FIG. 10: Convex hull and metastable phases of the B-Y-Zr ternary system. For binary structure types see Figs, I4I51HI Plotting 
symbols as in Fig. Energy units are meV/atom. 



portant energetic role for this system. Interestingly, this 
is the only structure among the stable and metastable B- 
Fe compounds (with the exception of cF116) in which B 
atoms do not have trigonal-prismatic environment. Each 
B atom has two other B atoms only 2.1A distant, forming 
ID chains along the c-direction. 

B atoms in the cF116 structure have the (0,5,4) 
Voronoi polyhedron; if we eliminate its smallest face (area 
0.35A 2 ), it converts to an (0,8,0) environment, in which 
each B atom has 8 Fe neighbors. This B environment 
is similar to that of BFe2 except there are no B-B near 
neighbors. The cF116 structure is also exceptional in an 
uneven distribution of the B atoms in Fe matrix: while 
some Fe atoms (sites 4a and 8c) have no B neighbors, 
site 32f has two B neighbors and 48h has three. 



environments; (ii) the (0,5,4), with 8 Fe and 1 B atom 
environments. These occur in roughly equal proportion. 
The trigonal prisms are consistent with the main struc- 
tures of the Fe-rich metastable phases, while the (0,5,4) 
environments are characteristic of the stable BFe2-tI12 
structure in which B-B neighbors occur. In fact we oc- 
casionally find a B atom with two B neighbors, resulting 
in local environments very close to BFC2. 

Diffraction data finds no B-B neighbors, but they 
occur robustly in our simulations and perhaps can serve 
as nucleation sites for crystalline BFe2- There is some 
controversy in the literature about the certainty with 
which B-B neighbors can be ruled out [5l| . 



3. Boron environments in Fe-rich ternaries 



2. Quenched samples of Fe$oB2o 

In the quenched samples, relaxed to the local minimum 
in cohesive energy at OK, the most common B environ- 
ments are: (i) trigonal prisms (0,3,6) with 9 Fe atom 



For B-Fe-Y, in the stable compounds BCeCo4.hP12 
and BFei 4 Y 2 .tP68, as well as B 2 Nd 3 Nii 3 .hP18 struc- 
ture (which in our calculation is unstable by just 
4meV/atom), the B environment is a trigonal prism with 
Fe at vertices and Y capping the rectangular faces. In- 



14 



terestingly, in the metastable B3Fe23Nd2.cI224 (this lies 
just 12 meV/atom above the tieplane) the trigonal prisms 
come in pairs, sharing rectangular faces and creating one 
B-B bond per pair. In contrast, in the cF116 structure, 
which is nearly stable in the B-Fe-Zr system, B atoms 
are surrounded by Fe atoms only. 



B. Structure of Fei7Y2 

The structure of the compound Fei7Y 2 is not pre- 
cisely known. Multiple structural variants have been 
observed, and the best structure refinements contain 
many partially occupied sites. The structures have close 
structural relationship to FesY.hP6 (CaCus proto- 
type), in which Y atoms center hexagonal columns of 
Fe atoms. The columns of Y atoms (c = 4.1 A) form 
a triangular lattice with edge length a = 4.9A. Start- 
ing from this structure, the Fei7Y2 family may be de- 
rived by (i) taking the superstructure defined by the 
vectors (1,-1,0) x (1,2,0) x (0,0,2); (ii) applying the 
substitutional rule Y— »2Fe. Neighboring Y atoms (sepa- 
rated by either a or c distances) should never be substi- 
tuted simultaneously. This rule enforces planar hexag- 
onal lattices of Y atoms, with apparent stacking de- 
grees of freedom [44|. The Th 2 Zni7.hR19 prototype (also 
known as a) takes the ABC stacking sequence (we de- 
note 4A bilayers by capital letters) with c a ~ 12A. The 
Nii7Tli2.hP38 prototype (also known as (3) takes the AB 
sequence with ~ 8A. The two reported refinements of 
Fei 7 Y 2 (hP44 ; 52j and hP80 H3) are apparently disor- 
dered versions of /3-Nii7Th 2 . 

Our calculation confirms small energy differences be- 
tween the stacking variants: we find AB = AC is 0.9 
meV/atom higher in energy than ABC, which in turn 
is 0.9 meV/atom higher in energy than the the ABAC 
sequence with ~16A stacking period. The ABAC se- 
quence, which is the best we have found, leads to frac- 
tional occupancy of some Wyckoff positions in qualitative 
agreement with the hP44 refinement |52(. 

Experimentally reported fractional occupancies, along 
with occupancies of the stacking sequences we studied, 
are reported in Table ITTT1 The Y^2Fe substitution rule 
together with our assumption of a disordered "A * A*" 
stacking sequence, constrains site occupancies p so that 
p(Yl)+p(Y2)=l, and p(Fel)+p(Fe2)=l. Thus the hP44 
refinement implies some Y vacancies, while the hP80 re- 
finement places extra Fe atoms at Fel+Fe2. We consid- 
ered these possibilities, but found both of them energet- 
ically unfavorable. Therefore we believe that Fei7Y2 is 
the correct stoichiometry, and the mismatches in occu- 
pancy factors are artifacts of the refinement, arising from 
stacking disorder. 

The absence of well-ordered crystalline samples further 
supports our proposal that the Fei7Y2 phase could be 
unstable at low temperatures. 



TABLE III: Fractionally occupied Wyckoff sites in Fei7Y2 
structures. First column labels sites as in the hP44 refine- 
ment |H3l . second column (/i) gives the number of equiva- 
lent atoms per unit cell. Stackings AB and AC are differ- 
ent crystallographic settings, but otherwise completely equiv- 
alent. Fe a sites appear when we register the hR19 structure 
of Th 2 Zni 7 (ABC stacking) with the hP38, hP44 and hP80 
structures. This atom is not present in either hP44 or hP80 
refinements. The final row reports AE in units of meV/atom. 



site 


/' 


AB 


AC ABAC ABC hP44 hP80 


Yl 


2 


1 


1/2 


3/4 


2/3 


0.41 


0.71 


Y2 


2 





1/2 


1/4 


2/3 


0.35 


0.12 


Y3 


2 


1 


1 


1 


2/3 


1.0 


1.0 


Fel 


4 





1/2 


1/4 


1/3 


0.28 


0.29 


Fe2 


4 


1 


1/2 


3/4 


1/3 


0.71 


0.86 


Fe a 


4 











1/3 








AE 


21.4 


21.4 


20.7 


21.0 







C. Glass formation 

By inspection of our cohesive energy data we can iden- 
tify the main crystal phases that are likely to compete 
with formation of the amorphous solid. In the vicinity of 
Fe8oB2o, the structure C6Cr23.cF116 can crystallize with 
almost no composition shift. However, it may be diffi- 
cult to nucleate and grow such a complex crystal type 
during a quench, so a simpler nearby structure such as 
Fe3B.oP16 may be favored. In the limit of slow cooling, 
phase separation into pure Fe and Fe2B will occur. In 
fact, all these structures are reported in annealed sam- 
ples of Fe-B glasses. 

Alloying with Zr can be advantageous because 
the large and strongly interacting Zr atoms diffuse 
slowly p54|. Due to its large size, Zr strongly destabilizes 
the Fe3B.oP16 structure. However, we see in table HVl 
that Zr actually tends to stabilize slightly the CgC^ 
structure, and also risks formation of the BFei4Nd2.tP68 
structure. Stabilization of the CgCr23 structure may be 
counteracted with replacement of Y for Zr. 

Another danger of alloying with Zr is formation of the 
highly stable binary B2Zr.hP3. Indeed, in the B-Fe-Zr 
ternary [l(j even Fe-rich liquids coexist with solid B2Zr. 
Choosing compositions with B content below the eutec- 
tic (17% B) can avoid B2Zr formation. Also, alloying 
with the less strongly interacting (but still very large) 
element Y can counteract this. However, it is not advan- 
tageous to alloy only with Y because 1) Y atoms diffuse 
more quickly than Zr as a result of their weaker binding; 
2) There are several Fe-rich Fe-Y binary structures (e.g. 
Fe5Y, Fei7Y 2 ) whose formation should be avoided. A 
reasonable composition that balances these difficulties is 
B 15 Fe75Y 3 Zr7. 

Recent experimental studies [5^, |5(j of Y addition to 
Fe-Zr-B-based glasses confirm its beneficial effects, while 
attributing them to different causes [5o| . 
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TABLE IV: Enthalpies of crystal phases competing with glass 
formation. 





Fe 3 B.oP16 


C 6 Cr 2 3.cF116 


BFei 4 Nd 2 .tP68 


quench 


B-Fe 

B-Fe-Y 

B-Fe-Zr 


-211 (17) 
-27 (253) 
-218 (81) 


-168 (18) 
-163 (74) 
-253 (16) 


+215 (268) 
-115 (0) 
-116 (50) 


-29 (151) 
-59 (186) 
-79 (221) 



structures but larger than the 4.6-4.7A 3 occurring in the 
metastable crystalline phases. Because the energies and 
geometry of the CCT models are close to the amorphous 
structure, we call the CCT models "amorphous approxi- 
mants" . 



D. Canonical-cell models of amorphous 
approximants 

Canonical cell tilings [57| (CCT) form networks of 
icosahedral cluster centers for models of icosahedral qua- 
sicrystals. They were motivated by the cubic "1/1" 
approximants of the quasicrystals, in which icosahedral 
clusters are located at the vertices of BCC lattice, and 
connected by 2-fold ("6") and 3-fold ("c") inter-cluster 
linkages. Linkages of 6-type are longer than c-type link- 
ages by the factor 2/-\/3=1.15. The be- network with 
global icosahedral symmetry and maximal density of 
clusters is a tiling of four canonical cells, A, B, C and D. 
An A-cell is the two-fold symmetric BCC-tetrahedron, B- 
cell is a skewed rectangular pyramid, C cell is a three-fold 
symmetric tetrahedron and D-ce\\ is a trigonal prism. 
The trigonal prism has equilateral triangular bases of b- 
type linkages and rectangular faces of b- and c-type link- 
ages. 

Consider a model amorphous structure in which the 
icosahedral clusters are replaced by "large" atoms, and 
the largest of the cells, the trigonal prism D, is deco- 
rated by a small atom in the center. Such a model ap- 
pears to be entirely plausible for B-Fe, since: (i) ideal B 
atom environments in the D-ce\\ produce very nearly the 
same ratio of B-Fe/Fe-Fe nearest neighbor distances as 
the BFe3 compounds; (ii) provided the density of Z?-cells 
is fixed by B atom content, the canonical cells models 
should yield optimal Fe-Fe connectivity via the b and c 
linkages, with b/c length ratio fixed at the BCC structure 
value 2 / -\/3. In the limit of large unit cell size approxi- 
mating an icosahedral quasicrystal, the B content of this 
model ranges from xb ~ 0.11 to 0.22, covering the best 
glass-forming composition range. 

As a convenient and simple example of the CCT model, 
we consider smallest cubic approximant model containing 
all kinds of cells, so called "3/2" tiling with 32 CCT nodes 
per cubic cell, and Pa3 space group. It contains 8 D-cells 
centered by B atoms (also 72 A-cells, and 32 each of B- 
and C-cells), and has composition Fe4B. This ideal model 
remains practically undistorted upon relaxation, and re- 
laxes to an energy ^210 meV/atom above the tie-line. 
B-Fe and Fe-Fe bond lengths in the model are similar to 
those found in the metastable crystalline compounds of 
similar composition. Performing Voronoi analysis we find 
that the Fe atom polyhedron volume is similar to that of 
the crystals, but the B atom polyhedron volume is larger. 
In the CCT model the B atom polyhedron volume is 5.03 
A 3 , similar to the volume of 4. 9 A 3 in relaxed amorphous 



V. CONCLUSIONS 

In conclusion, we present a method for the calcula- 
tion of low temperature (T=0K) alloy phase diagrams 
and apply it to the study of the B-Fe-Y-Zr quaternary 
system. A key distinguishing feature of our approach is 
the establishment of two databases: one from which we 
draw promising structures observed in similar chemical 
systems; the other in which we record our calculated co- 
hesive energies which can be quickly converted into en- 
thalpies of formation. The first database allows us to 
propose and evaluate candidate structures even in alloy 
systems that have not been previously studied. The sec- 
ond allows us to quickly add a new chemical element, and 
re-use, for example, all our B-Fe-Zr data in the study of 
the B-Fe-Y-Zr quaternary. 

As a result we have examined two new ternary sys- 
tems (B-Y-Zr and Fe-Y-Zr) and the B-Fe-Y-Zr quater- 
nary. We found certain binary phases extending into the 
new ternary systems, but no ternary phases extending 
into the quaternary. So far we have not discovered any 
stable quaternary structure. Even in previously studied 
binary and ternary systems we find some new results in- 
cluding proposed structures for previously unsolved com- 
pounds. 

The broad agreement between our calculations and 
experimentally reported phase diagrams demands that 
special attention be paid where disagreements exist. 
These disagreements fall into certain categories: Un- 
certain reported structure found to have high energy 
(e.g. BZr.cF8); Well established experimental structure 
found to have high energy (e.g. Fei7Y 2 ); Structure cal- 
culated to be stable not present in published diagram 
(e.g. B 4 CoZr in B-Fe-Zr and Co 7 Y 2 in Fe-Y). These dis- 
agreements warrant further study, both theoretical and 
experimental. 

On the subject of glass formation, the main motivation 
for this study, we identify important crystalline competi- 
tors to glass-formation and illustrate how they can be 
destabilized by the addition of appropriately chosen large 
atoms. 
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